Random matrix analysis of complex networks 
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We study complex networks under random matrix theory (RMT) framework. Using nearest- 
neighbor and next-nearest-neighbor spacing distributions we analyze the eigenvalues of adjacency 
matrix of various model networks, namely, random, scale-free and small-world networks. These 
distributions follow Gaussian orthogonal ensemble statistic of RMT. To probe long-range correlations 
in the eigenvalues we study spectral rigidity via A3 statistic of RMT as well. It follows RMT 
prediction of linear behavior in semi-logarithmic scale with slope being ~ 1/tt 2 . Random and 
scale-free networks follow RMT prediction for very large scale. Small-world network follows it for 
sufficiently large scale, but much less than the random and scale-free networks. 
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I. INTRODUCTION 

Random matrix theory (RMT), initially proposed to 
explain statistical properties of nuclear spectra, had suc- 
cessful predictions for the spectral properties of different 
complex systems such as disordered systems, quantum 
chaotic systems, large complex atoms, etc., followed by 
numerical and experimental verifications in the last few 
decades 0- Quantum graphs, which model the sys- 
tems of interest in quantum chemistry, solid state physics 
and transmission of waves, have also been studied under 
the RMT framework 0] ■ Recently, RMT has been shown 
to be useful also in understanding the statistical proper- 
ties of empirical cross-correlation matrices appearing in 
the study of multivariate time series of followings: price 
fluctuations in stock market [1, [B| , EEG data of brain , 
variation of different atmospheric parameters 0], etc. 

In our previous studies 0, 0] complex networks have 
been analyzed under RMT framework. These works con- 
sider nearest-neighbor spacing distribution (NNSD) of 
eigenvalues spectra of adjacency and Laplacian matri- 
ces of various extensively studied networks. The NNSD 
gives probability for finding neighboring eigenvalues with 
a given spacing, and it follows two universal properties 
depending upon underlying correlations among the eigen- 
values. For the correlated eigenvalues, NNSD follows 
Gaussian orthogonal ensemble (GOE) statistics of RMT, 
whereas it follows Poissonian statistics for the uncorre- 
cted eigenvalues. One of the main advantages of RMT 
approach is that depending on the nature of eigenval- 
ues correlations one can separate system dependent part 
from random universal part, which are intermingled due 
to the complexity of the system @, 0, la, la 01 • RMT 
analysis for the various networks shows that the NNSD 
of complex networks also follow universal GOE statistics 
of RMT [8] . This finding suggests that different results of 
GOE statistics, which have successfully been applied to 
understand the systems coming from various fields start- 
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ing from nuclei to the stock-market, can be applied to 
study networks as well. Our earlier works concentrate on 
the NNSD studies of networks. NNSD carries informa- 
tion for the correlation between two adjacent eigenvalues, 
but do not tell about the correlation between two far- 
off eigenvalues. Therefore, even though NNSD follows 
GOE statistics of RMT, other properties may show de- 
viations, which suggests that one can not rely on NNSD 
results exclusively. To probe for long-range correlations 
as well, current work considers spectral rigidity test via 
well known A 3 -statistic of RMT. It is found that the spec- 
tral rigidity of the complex networks follows RMT pre- 
diction, with scale depending upon the properties of the 
networks. Present work also analyzes the neii-nearest- 
neighbor spacing distribution (nNNSD) of the adjacency 
matrix of the networks. 

The paper is organized as follows: following the above 
introduction, Sec. HT1 explains various aspects of complex 
networks studies. Sec. IIIII describes some basics of RMT 
relevant to our studies. Sec. HVI illustrates the RMT anal- 
ysis for various model networks, namely; random, scale- 
free and small world. The NNSD is the most widely stud- 
ied property in random matrix literature, therefore, this 
section includes NNSD results for the above mentioned 
model networks 0], and presents results for nNNSD and 
the A3-statistic of these networks. Finally, Sec. [V] dis- 
cusses and summarizes results with some possible future 
directions. 



II. COMPLEX NETWORKS 

Last 10 years have witnessed a rapid advancement in 
the studies of complex networks. The main concept of 
the network theory is to define complex systems in terms 
of networks of many interacting units. Few examples 
of such systems are interacting molecules in living cell, 
nerve cells in brain, computers in Internet communica- 
tion, social networks of interacting people, airport net- 
works with flight connections, etc (lol. Illl H3 |. In the 
graph theoretical terminology, units are called nodes and 
interactions are called edges fl3| . Various model net- 
works have been introduced to study the behavior of 
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complex systems having underlying network structures. 
These model networks are based on simple principles, still 
they capture essential features of the underlying systems. 



A. Structural properties 

In random graph model of Erdos and Renyi any two 
nodes are randomly connected with probability p [3]. 
This model assumes that interactions between nodes are 
random. Recently, with the availability of large maps of 
real world networks, it has been observed that the ran- 
dom graph model is not appropriate for studying the be- 
havior of real world networks. Hence many new models 
have been introduced. Watts and Strogatz proposed a 
model, popularly known as 'small- world network', which 
has properties of small diameter and high clustering [l5[ . 
Moreover, this model network is very sparse : network 
with a very few number of edges, another property shown 
by many real-world networks. In addition to the above 
mentioned properties, Barabasi and Albert show that 
degree distributions of many real-world networks have 
power-law. This implies that some nodes are much more 
connected than the others [lfij]. Barabasi- Albert's scale- 
free model and Watts-Strogatz's small-world model have 
contributed immensely in understanding evolution and 
behavior of the real systems having network structures. 
Following these two new models came an outbreak in the 
field of networks. These studies show that real world 
networks have coexistence of randomness and regularity 

[a nam. 



B. Spectral properties 

Apart from the above mentioned investigations which 
focus on direct measurements of the structural properties 
of networks, there have been lot of studies demonstrat- 
ing that properties of networks or graphs could be well 
characterized by the spectrum of associated adjacency 
matrix A [19J. For an unweighted graph, it is defined in 
the following way : Aij = 1, if i and j nodes are con- 
nected and zero otherwise. For an undirected network, 
this matrix is symmetric and consequently has real eigen- 
values. Eigenvalues give information about some basic 
topological properties of the underlying network [n| [2(| • 
Spectral properties of networks have also been used to 
understand some of the dynamical properties of interact- 
ing chaotic units on networks, for example largest eigen- 
value of adjacency matrix determines transition to the 
synchronized state 21]. The distribution of the eigen- 
values of a matrix having finite probability p of nonzero 
Gaussian distributed random elements per row, follows 
Wigner semicircular law in the limit p — > 1. For very 
small p also, which corresponds to the sparse random 
matrix, one gets semicircular law with several peaks at 
different eigenvalues (22| . 



With the increasing availability of large maps of real- 
world networks, analyses of spectral densities of adja- 
cency matrix of real- world networks and model networks 
having real- world properties have also begun [23|, [U, HH, 
Hi] . These analyses show that the matrix constructed 
by zero and one elements corresponding to a unweighted 
random network, also follows Wigner semicircular law 
[23| with degeneracy at A = 0. Small- world model net- 
works show very complex spectral density with many 
sharp peaks [26| , while the spectral density of scale-free 
model networks exhibits so called triangular distribution 
[H HS HI]. Spectral density and NNSD of the random 
matrices constructed by zero and one elements have been 
studied extensively in the Ref. [27| . These studies show 
that NNSD of the random matrices follow GOE distri- 
bution of RMT. 



III. RANDOM MATRIX STATISTICS 

In the random matrix studies of eigenvalues spectra, 
one has to consider two kinds of properties : (1) global 
properties, like spectral density or distribution of eigen- 
values p(A), and (2) local properties, like eigenvalue fluc- 
tuations around p(A). Among these, the eigenvalue fluc- 
tuations is the most popular one. This is generally ob- 
tained from the NNSD of the eigenvalues. The eigenval- 
ues of network are denoted by Ai, i = 1, . . . , N, where 
N is the size of the network and A^ < Aj+i, Vi. In order 
to get universal properties of the eigenvalue fluctuations, 
one has to remove the spurious effects due to the varia- 
tions of spectral density and to work at constant spectral 
density on the average. Thereby, it is customary in RMT 
to unfold the eigenvalues by a transformation Ai = -/V(A;), 

where N(X) — J*. p(A') dX' is the averaged integrated 

eigenvalue density Since analytical form for N is 
not known, we numerically unfold the spectrum by poly- 
nomial curve fitting. Using the unfolded spectrum, we 
calculate the nearest-neighbor spacings as 

s l ' — A,;+i — A^ 

and due to the above unfolding, the average nearest- 
neighbor spacings (si) becomes unity, being independent 
of the system. The NNSD P(si) is defined as the proba- 
bility distribution of these s^' 's. In case of Poisson statis- 
tics, 

P(si) = exp(-ai); (1) 

whereas for GOE 

P( Sl ) = J Sl exp(-^r). ( 2 ) 

For the intermediate cases, NNSD is described by Brody 
formula [H]: 

P ( Sl ) = A^exp(-B^ +1 ) , (3a) 
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where A and B are determined by the parameter f3 as 
follow : 



A = (1 + /3)aanda 



/3 + 2 
/3 + 1 



0+1 



(3b) 



This is a semiempirical formula characterized by the sin- 
gle parameter (3, popularly known as Brody parameter. 
(3 = 1 corresponds to the GOE statistics and (3 = cor- 
responds to the Poisson statistics. 

Apart from NNSD, the nerrf-nearest-neighbor spacings 
distribution (nNNSD) is also used to characterize the 
statistics of eigenvalues fluctuations. We calculate this 
distribution P(s2) of next-nearest-neighbor spacings 



,(0 



(\+2 — Ai)/2 



(4) 



between the unfolded eigenvalues. Factor of two at the 
denominator is inserted to make the average of next- 
nearest- neighbor spacings (S2) unity. According to Ref. 
1] , the nNNSD of GOE matrices is identical to the NNSD 
of Gaussian symplectic ensemble (GSE) matrices, i.e., 



P(S2) 



3 6 7T 3 



s\ exp 



64 
9^ 



(5) 



The NNSD and nNNSD reflect only local correlations 
among the eigenvalues. The spectral rigidity, measured 
by the A 3 -statistic of RMT, gives information about the 
long-range correlations among eigenvalues and is more 
sensitive test for RMT properties of the matrix under 
investigation [l|, |29|. In the following we describe the 
procedure to calculate this quantity. 

The A3-statistic measures the least-square deviation of 
the spectral staircase function representing the averaged 
integrated eigenvalue density N(X) from the best straight 
line fitting for a finite interval L of the spectrum, i.e., 



A 3 (L; x) = — min 

L/ a,b 



[N{\) - aX 



b] 2 dX 



(6) 



where a and b are obtained from a least-square fit. Av- 
erage over several choices of x gives the spectral rigidity 
A 3 (L). For the Poisson case, when the eigenvalues are 
uncorrelated, A 3 (L) = L/15, reflecting strong fluctua- 
tions around the spectral density p(X). On the other 
hand, for the GOE case, A 3 (Z) depends logarithmically 



on L, i.e. 



A 3 (L) 



InL. 



(7) 



IV. RESULTS 

In the following we present results for the ensemble av- 
eraged NNSD, nNNSD and A 3 statistic of random, scale- 
free and small- world networks. 
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FIG. 1: (Color online) Nearest-neighbor spacings distribution 
(NNSD) -P(si) of the adjacency matrices of different networks 
[(a) random network, (b) scale-free network, and (c) small- 
world network]. All follow GOE statistics. The histograms 
are numerical results and the solid lines represent fitted Brody 
distribution (Eq.(3)). All networks have TV = 2000 nodes and 
an average degree k — 20 per node. Figures are plotted for 
average over 10 random realizations of the networks. Insets 
show respective spectral densities p(A). 



A. Random network 

First we consider an ensemble of random networks gen- 
erated by using Erdos-Renyi algorithm. Starting with 
N = 2000 nodes random connections between pairs of 
nodes are made with probability p. The average de- 
gree of the graph is k = 2n/N = p(N — 1) ~ pN. 
There exists a critical probability p c (N) for which one 
gets a large connected component. The degree distri- 
bution of this random graph is a binomial distribution 
P(k) = C k N _ lP k {l - p) N - 1 ~ k . For p = 0.01, this 
method yields a connected network with average degree 
p x N = 20. Note that for very small value of p one gets 
several disconnected components. In this study choice of 
p is high enough to give a connected component typically 
spanning all the nodes. 

We calculate the eigenvalues spectrum of network gen- 
erated according to the above algorithm. First the eigen- 
values are unfolded by using the technique described in 
Sec. IIIII This method yields the eigenvalues with con- 
stant spectral density on the average. These unfolded 
eigenvalues are used to calculate NNSD. The same pro- 
cedure is repeated for an ensemble of the networks gen- 
erated for different random realizations. Note that p is 
always chosen such that algorithm generates a network 
with average degree k = 20. FigQJa) plots ensemble av- 
erage of NNSD. By fitting this ensemble averaged NNSD 
with the Brody formula given in Eq.Q we get an esti- 
mation of the Brody parameter j3 = 0.9786 ~ 1. This 
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FIG. 2: (Color online) Next-nearest-neighbor spacings distri- 
bution (nNNSD) P(s2) of the adjacency matrices of different 
networks [(a) random network, (b) scale-free network, and (c) 
small-world network] is compared with the nearest-neighbor 
spacings distribution (NNSD) of GSE matrices. Figures are 
plotted for average over 10 realization of the networks. All 
networks have N = 2000 nodes and an average degree k = 20 
per node. 



FIG. 3: (Color online) Aa(L) statistic for eigenvalues spectra 
of the random network. The circles are numerical results and 
the solid curve is GOE prediction of RMT. Inset shows the 
As(L) in semi-logarithmic scale, in this scale it has the slope 
0.0978. Figure is plotted for average over 10 realizations of 
the networks. All networks have N — 2000 nodes and an 
average degree k = 20 per node. 



value of Brody parameter clearly indicates the GOE be- 
havior of the NNSD [Eq. ©]. Inset of Fig. QJa) shows 
corresponding spectral density which follows well known 
Wigner's semicircular distribution. The same unfolded 
eigenvalues are used to calculate nNNSD. For this we cal- 
culate next nearest neighbor spacings as given in Eq. (U|) 
and plot their distribution in Fig. [UJa) . It can be seen 
from the figure that the nNNSD agrees well with the 
NNSD of GSE matrices as given in Eq. ©. 

As explained in the introduction that NNSD and 
nNNSD only tell about the short range correlations 
among the eigenvalues. Therefore, to probe for the long 
range correlations we study A^(L) statistic of the spec- 
trum of this network. A^(L) is calculated following 
Eq. ©. Fig. [3] plots this statistics for the same ensemble 
as used for the NNSD and nNNSD calculations above. It 
can be seen that A3 (L) statistic agrees very well with the 
RMT prediction, given by Eq. 0, upto very large value 
of L, i.e., L ~ 300. Inset of this figure shows the same 
in semi-logarithmic scale. Here one can see the expected 
linear behavior of A^(L) with slope of 0.0978 which is 
very close to the RMT predicted value l/?r 2 ~ 0.1013 
[Eq.©]. 

Note that here an ensemble of ten networks of dimen- 
sion N = 2000 is considered. Statistical properties of 
eigenvalues spectra of members of this ensemble have 
very small deviations from each other and hence justi- 
fying ensemble averaging calculations (30j . Each indi- 
vidual network in the ensemble follows random matrix 
predictions with very good accuracy, however to make 
the statistical analysis more credible, we present the re- 



sults for an ensemble of ten networks. Here we would 
like to mention that an ensemble of networks of much 
smaller dimensions, say N — 100, has been studied as 
well and it follows GOE predictions of RMT. However, 
for this case, many more realizations are required to get 
good accuracy. 



B. Scale-free network 

Scale-free network is g enerated by using the model of 
Barabasi and Albert [161 ]. Starting with a small number, 
mo of the nodes, a new node with m < mo connections 
is added at each time step. This new node connects with 
an already existing node i with probability 7r(fc,-) cx fcj, 
where fcj is the degree of the node i. After r time steps the 
model leads to a network with N = t + tuq nodes and mr 
connections. This model generates a scale-free network, 
i.e., the probability P(k), that a node has degree k decays 
as a power law P(k) ~ k~ x , where A is a constant and for 
the type of probability law ir(k) used here A = 3. Other 
forms for the probability 7r(fc) are also possible which give 
different values of A. However, the results reported here 
are independent of the value of A. 

Using the above algorithm an ensemble of scale-free 
networks of size N — 2000 and average degree k — 20 is 
generated. To calculate NNSD, nNNSD and A 3 for the 
spectra of this ensemble, we follow the same procedure 
as described in the previous section. Fig. []Jb) shows 
that the NNSD of scale-free network follows GOE with 
(3 = 0.9626 ~ 1. Inset of this figure plots the spectral 
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FIG. 4: (Color online) Aa(L) statistic for eigenvalues spectra 
of the scale-free network. The circles are numerical results and 
the solid curve is the GOE prediction of RMT. Inset plots the 
As(L) in semi-logarithmic scale, in this scale it has the slope 
0.0975. Figure is plotted for the average over 10 realizations 
of the networks. All networks have N = 2000 nodes and an 
average degree k — 20 per node. 



density of this network. Fig.[5Jb) shows that the nNNSD 
of the adjacency matrix of this network agrees well with 
the NNSD of the GSE matrices. Fig. H shows the A 3 (L) 
statistic for the adjacency matrix of scale-free network. 
Here we see that the As(L) statistic for the scale- free 
network agrees very well with the RMT prediction for 
very large L, i.e., L ~ 150, and deviations are seen only 
after L = 150. Inset of this figure shows the expected 
linear behavior of A^(L) in semi-logarithmic scale for L < 
150 with the slope of 0.0975, a value very close to the 
RMT predicted value 1 /it 2 . 

Universality of NNSD and nNNSD for random and 
scale-free networks seems to give the impression that 
these networks have same amount of randomness but A3 
results tell that the scale-free network is not as much 
random as the random network. This is obvious from 
their construction algorithms as well, but A3 statistics is 
capturing this property which is a very important result. 
The finding also suggests that scale-free networks have 
some specific features that cannot be modeled by RMT. 
It may be noted that one can generate scale- free networks 
by using other algorithms as well [3jJ, |32| , for these net- 
works also spacing distributions and spectral rigidity re- 
sults will have qualitatively similar behaviors, except that 
the range of agreement of L with the random matrix pre- 
diction would depend upon the amount of randomness in 
the networks. 



C. Small- world network 

Small-world networks are constructed using the follow- 
ing algorithm of Watts and Strogatz [l5[ . Starting with 
an one-dimensional ring lattice of N nodes in which every 
node is connected to its k/2 nearest neighbors, we ran- 
domly rewire each connection of the lattice with the prob- 
ability p such that self-connections and multiple connec- 
tions are excluded. Thus p = gives a regular network 
and p — 1 gives a completely random network. The typi- 
cal small- world behavior is observed around p — 0.005 [8]. 
For N = 2000 and average degree k = 20, an ensemble of 
ten different realizations of the network are generated. 

Again the same procedure as described in Sec. IIV Al 
has been used to calculate NNSD, nNNSD and A 3 for the 
spectra. Figfljc) shows that the NNSD of this network 
again follows GOE statistics with j3 very close to 1, i.e., 
(3 = 1.0136. The inset shows that the corresponding spec- 
tral density is complicated with several peaks. One peak 
is always at A = 0. The exact positions of other peaks 
may vary but overall form of spectral density remains 
similar. Fig. [^c) plots the nNNSD of adjacency matrix 
of small-world network. It can be seen that the nNNSD 
agrees well with the NNSD of GSE matrices. Fig.[5]shows 
the As(L) statistic for the spectrum of adjacency matrix 
corresponding to the small-world network with p — 0.005. 
Inset of this figure shows the expected linear behavior of 
A^(L) in semi- logarithmic scale for L < 30 with slope of 
0.1024, a value very close to the RMT predicted value 
1/ir 2 . It can be seen here that A 3 (i) statistics for the 
small- world network agrees very well with the RMT pre- 
diction for sufficiently large L, i.e., L ~ 30, but much 
less than the same for random and scale-free networks, 
implying that besides randomness, small- world network 
has specific features also. This again suggests that the 
behavior of A3 statistics can be used to understand the 
amount of randomness in the networks. More specifically 
deviation from the GOE predicted behavior corresponds 
to the system specific features in the network. 

Note that in this paper, results for networks with the 
average degree 20 are presented. We have studied sparser 
((k) < 20) and denser networks ((k) > 20 to (k) ~ N) 
as well. Same universal behavior are found for these net- 
works as far as there exists a certain amount of random- 
ness, i.e., presence of some minimal random connections 
among the nodes. There exists problems with very sparse 
networks as of average degree two and very dense net- 
works as of degree ~ N. For sparse networks ((fc) ~ 2) 
sometimes one can get several degeneracies in eigenvalues 
[22| . In this case, one has to first get rid of the degen- 
eracies to conclude anything under the RMT framework 
& 

Similarly, for dense networks, universal spacing distri- 
bution are observed till very large value of average de- 
gree. For (k) ~ N, largest eigenvalue has very high value 
than the rest of N — 1 eigenvalues which are very close 
to each other, becoming equal in the limiting case of all 
to all connections (k) = N — 1. For example, random 
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FIG. 5: (Color online) Aa(L) statistic for the eigenvalues 
spectra of the small- world network. The circles are numer- 
ical results and the solid curve is GOE prediction of RMT. 
Inset plots the As(L) upto L — 30 in semi- logarithmic scale, 
in this scale Aa(L) has the slope 0.1024. Figure is plotted for 
the average over 10 realizations of the networks. All networks 
have N — 2000 nodes and an average degree k = 20 per node. 



networks with p = 0.95 (which means that the network 
has 95% of maximum possible connections), also follow 
RMT predictions of universal spacing distributions till 
very large scale. As the number of connections are in- 
creased further one starts getting degenerate eigenvalues 
and for p = 0.999 high degeneracies at various values 
(such as A = —1,0,1) are observed keeping it trivially 
out from the RMT studies. 



V. DISCUSSIONS 

We use RMT to study complex networks and show 
that in spite of spectral densities of the adjacency ma- 
trices being different for different networks, their eigen- 
value fluctuations are same and follow the GOE statistics 
of RMT. We attribute this universality to the existence 
of minimal amount of randomness in all these networks 
and show that randomness in the network connections 
can be quantified by the Brody parameter. In addition 
to the NNSD, we present the results of nNNSD and spec- 
tral rigidity via A 3 -statistic of RMT. The nNNSD of the 
eigenvalues of these model networks are identical to the 
NNSD of GSE matrices which again agrees with RMT 
prediction given in Eq. ([5]). NNSD and nNNSD sug- 
gest that there exists short range correlations among the 
eigenvalues. Spectral rigidity test shows that the A 3 (L) 
statistics follows random matrix predicted linear behav- 
ior in semi-logarithmic plot for sufficiently large scale L 
with slope being ~ 1/7T 2 [see, Eq.®], suggesting long- 
range correlations among the eigenvalues. Above find- 



ings show that statistics of the bulk of eigenvalues of the 
model networks is consistent with those of a real sym- 
metric random matrix and deviation from this could be 
understood as a system dependent part. 

Universal GOE behavior of NNSD and nNNSD tell 
that the networks are sufficiently random, or there ex- 
ists minimal amount of randomness required to introduce 
the correlations among the neighboring eigenvalues. The 
A3 analysis seems to characterize the level of random- 
ness in networks depending on the range of correlations 
among eigenvalues. A3 analysis of the random network 
follows RMT prediction for very long range of L, which 
is not very surprising as random network follows RMT 
at each level starting from semi-circular density distribu- 
tions. However interestingly scale- free and small- world 
networks also follow RMT for sufficiently large value of L. 
Beyond this value of L, deviation in the spectral rigidity 
is seen, indicating a possible breakdown of the universal- 
ity. This is quite understandable as small-world network 
is highly clustered and scale-free network also has spe- 
cific features like hubs, so it is natural that they are not 
as random as the random network. But it is interesting 
to realize that A3 statistics rightly captures this informa- 
tion. Moreover, small-world network is generated exactly 
at small-world transition by using Watts and Strogatz 
algorithm which yields network with very high cluster- 
ing coefficient and very less number of random connec- 
tions. Results presented in this paper show that these 
very small number of random connections make network 
sufficiently random to introduce the correlations among 
the eigenvalues for the sufficiently long range. 

According to the many recent studies, randomness in 
connections is one of the most important and desirable 
ingredients for the proper functionality or the efficient 
performance of systems having underlying network struc- 
tures. For instance, information processing in brain is 
considered to be because of random connections among 
different modular structures dU . We feel that the role of 
random connections, and behavior and evolution of such 
systems can be studied better under the RMT frame- 
work. Also this RMT approach may be used to detect 
the connections most responsible to increase the complex- 
ity of networks. For example effect of oxygen molecule 
on biochemical network of a metabolic system is recently 
studied and is shown to increase the complexity of the 
system leading to a major transition in the history of life 

In summary, we use RMT to analyze spectra of com- 
plex networks and show that these networks follow uni- 
versal GOE statistics. These results tell that random 
matrix theory, a very well developed branch of Physics, 
can be applied to the complex networks studies. So far 
we have only concentrated on the model networks stud- 
ied vastly in the recent literature, providing a basis to 
the random matrix analysis of networks. Future investi- 
gations would involve studies of real- world networks [3(| . 
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